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We numerically construct charged boson stars in asymptotically AdS spacetime. We find an 
intricate phase diagram dominated by solutions whose main matter contribution are alternately 
provided by the scalar field or by the gauge field. 
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Introduction 

The study of boson stars dates back to the late 1960's 
with the work of Kaup [l| (see also @) who, inspired 
by ideas of Wheeler, constructed Klein- Gordon geons. 
These objects and their generalizations have found nu- 
merous applications in and beyond general relativity. 
Some of the classical reviews in the subject include: [3rtZ| . 

In this letter we construct charged boson stars in 
asymptotically Anti-de-Sitter (AdS) spacetimes and in- 
vestigate their properties. Namely, we consider particle- 
like solutions of a complex scalar field coupled to gravity 
and a Maxwell field in the presence of a negative cos- 
mological constant. We have three main motivations to 
study these objects. 

First, it is intrinsically interesting to understand 
particle-like solutions in asymptotically AdS space times 
to enhance and test our intuition of highly symmetric so- 
lutions of Einstein gravity with a negative cosmological 
constant. For example, the role of boundary conditions 
necessary to define dynamics in AdS is one aspect that 
constantly challenges our Minkowski-based intuition. 

Our second motivation is also linked to understand- 
ing dynamics in AdS but it is more concretely related 
to gravitational collapse in AdS. It has recently been es- 
tablished that AdS is unstable under arbitrarily small 
perturbations with respect to black hole formation Q. 
This result has been discussed and elaborated upon for 
massless fields Given previous history with crit- 

ical collapse of massive 14] and Yang-Mills fields fl5j j . 
it is possible that the phase diagram of critical collapse 
in asymptotically AdS spacetimes gets modified by the 
introduction of mass and other fields. In particular, the 
existence of boson stars might prevent a direct channel 
to black hole formation. 

Thirdly, solutions of a scalar field in Einstein-Maxwell 
gravity in asymptotically AdS spacetimes have proven 
a fruitful ground for applications of the AdS/CFT cor- 
respondence to various situations in condensed matter 
physics (see, for example, lB- 18|). 

Charged boson stars 

We consider a massive charged complex scalar field, <j>, in- 
teracting with electromagnetism and minimally coupled 



to Einstein gravity with a negative cosmological constant, 
S = I d 4 xyj—g 



— - — (R - A) - -F^F^ 



-(fl^)(D"f)- m ^> 



(1) 



D, = 



^-iqAp, F flv = V l± A u -V u A lil , (2) 



where G is Newton's constant. Since we are interested in 
stationary, spherically symmetric solutions, we write the 
metric in Schwarzschild form 



ds z 



-e 2u dt 2 + e 2v dp 2 + p 2 dn 2 2 . 



(3) 



Here p is the radial coordinate, and we take u — u(p) and 
v = v(p). We then take the scalar to be time-harmonic 



1 



V(p), 



(4) 



where o~(p) is a real function of p. Since <f> is electri- 
cally charged, it will source the electromagnetic field, and 
hence we turn on a scalar potential, A t — A t (p). 

At this point, a few quick comments are in order. 
Firstly, the Maxwell equation arising from ([1]) takes the 
standard form V^i 7 ^ = qj u , where 



Ju 



(5) 



is the conserved particle number current. As a result, 
the total charge of the boson star is given by Q = qN 
where N is the conserved particle number. Secondly, 
the time-dependence of <f> may be removed by a gauge 
transformation of the form A t — > A t — oj/q along with 
(\> — > e~ lut (f), so only the combination ui — qAt is physical. 
This is in contrast with the standard (ungauged) boson 
stars, where ui has an intrinsic meaning. 

It is straightforward to derive the coupled equations of 
motion corresponding to the above spherically symmet- 
ric ansatz. When developing the numerical solutions, we 
take as input the mass m and charge q of the scalar field 
as well as the cosmological constant A, and scale by New- 
ton's constant G when appropriate in order to work with 
dimensionless quantities. We also introduce the gauge 
invariant combination 



A(p) = (uj - qA t (p))/(iO - qA t (0)). 



(6) 
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To define our boundary value problem, we first de- 
mand regularity near the origin, p = 0. This leads to the 
following conditions: 

u(Q) = u Q , u(0)=0, cr(0)=cr o , o-'(0) = 
1(0) = 1, l'(0) = 0. (7) 

At asymptotic infinity, p — ¥ oo, AdS boundary conditions 
imply that the matter fields behave as 

a( P ) = a lp - A + a 2 p A ~ 3 , A = \ - ^- + (mL)\ 
A t (p) = ao + mp- 1 . (8) 

Our task is to look for solutions with normalizable modes 
at infinity, namely o~ 2 = 0. Therefore we have a boundary 
value problem which we solve numerically using shooting 
techniques. We determine a 2 = with a precision of 
10~ 20 . In our minimizing algorithm we shoot by changing 
do of the initial data for a given uo; we use (mL) = 10 
throughout. 

Properties of the solutions 

In the absence of a cosmological constant, there is a crit- 
ical charge, g cr it, above which the star does not exist in 
flat space. In the Newtonian limit the value of q c tit is 
obtained by comparing the gravitational attraction with 
the electrostatic repulsion of an elementary particle in- 
teracting with the star; the result is Gm 2 = q 2 rit / (47r). In 
the context of general relativity this result is only slightly 
modified. However, for boson stars in AdS, the value of 
g cl -it above which a star does not exist is substantially 
larger. Moreover, the phase space of solutions is quite 
different. 

Intuitively we can understand the increase of g cr it in 
the presence of a negative cosmological constant as fol- 
lows. Once the electrostatic force has overcome the grav- 
itational pull, the only force standing in the way of 
charged particles flying away is the pull generated by 
the cosmological constant. Therefore for large enough 
charges the main mechanism supporting the star is the 
balance between the gauge forces and the cosmological 
constant. 

For small charge q, the boson stars in AdS resemble 
those in flat space. The regular zero-node solutions have 
a smooth scalar profile, with maximum particle density 
at the core, and a gradual fall off as a function of p. 
Since the scalar field dominates the energy density, we 
denote these regular solutions as scalar dominated. For 
sufficiently large charge q, on the other hand, we find a 
new type of solution where the gauge field contribution 
dominates the total energy of the system. We naturally 
denote these as gauge field dominated. 

An example of the regular and new solutions is given 
in Fig. [TJ The regular scalar dominated solution is shown 
on the left, while the new gauge field dominated solution 
is on the right. There are three key properties that dis- 
tinguish the new solution from the regular one: (i) The 
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FIG. 1: The scalar field profile ff(p) (top) and energy con- 
tributions (bottom) for regular and new solutions. The left 
panel (q = 0.56) is for a standard charged boson star; the 
right panel (q — 0.6) is for a new solution. 



new solution has a sharper surface, as defined by the pro- 
file of the scalar field; (ii) The gauge field is concentrated 
near the surface of the star, instead of distributed in the 
interior as in the standard case; (Hi) The contribution 
to the total mass of the star is dominated by the gauge 
field. 

We now consider the mass of the boson star as a func- 
tion of its core density ct(0) and as a function of the 
particle number N. Since the interplay between regular 
and new solutions depends on the charge q, we define 
two critical charges, q\ and q 2 . For q < qi only the 
scalar dominated solution is possible. For qi < q < q 2 , 
we enter a transition region where both the scalar dom- 
inated and gauge field dominated solutions exist as dis- 
tinct branches. Finally, for q > q 2 , the scalar and gauge 
field dominated solutions merge. 

The mass as a function of c(0) and as a function of 
N for an intermediate value of the charge, q\ < q < q 2 , 
is shown in Fig. [2] The green curve represents the regu- 
lar scalar dominated solution, while the red curve cor- 
responds to the gauge field dominated solution. We 
have also shown the one-node solution as the blue curve. 
While this is ordinarily discarded as being an excited 
state, here we wish to note its proximity to the gauge 
field dominated solution. For each of these solutions, the 
right panel shows three curves of the same color. These 
three curves correspond to the total mass and the sepa- 
rate scalar and gauge field contributions to the mass. For 
the scalar dominated (green) curves, the scalar contribu- 
tion to the mass is the larger one, while for the gauge 
field dominated (red) curves, the gauge field contribu- 
tion to the mass is larger. Note that there is a gap on the 
right panel of Fig. [3] This suggests some sort of transi- 
tion between scalar dominance on the left and gauge field 
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dominance on the right. 
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FIG. 2: Mass as a function of scalar field at origin (left) and 
as a function of particle number (right) in the intermediate q 
regime. 

When the charge q is increased beyond the critical 
value qi we observe a very different behavior of the mass 
as a function of c(0) and as a function of N, as shown in 
Fig- El There is no longer a gap in the right panel and we 
see clearly the change in the contribution to the total en- 
ergy coming from the scalar-dominated region prevalent 
at small values of N and the gauge-dominated region at 
large values of N. The blue curve again corresponds to a 
one-node solution. This figure should be read as a merg- 
ing of the previous one in the cases where the charge has 
increased, that is, increasing the charge narrows the gap 
between the two kind of solutions which is evident in the 
region of qi < q < q2- 

The three green curves on the right panel of Fig. [3] arc 
total energy, the contribution from the scalar field and 
the contribution from the gauge field respectively. This 
graph should be read as a merging of branches: now the 
green curve is a mixture of the normal branch for low N 
and gauge-dominated branch for large N; the red curves 
are only in the areas of overlap in the previous graph. 
The most prominent feature of Fig. [3] is the crossing of 
the scalar and gauge contributions to the total energy. 

In the left panels of Figs. [2] and [3] we have also plotted 
a branch corresponding to solutions with one node in the 
amplitude of the scalar field profile, cr(p); these are tradi- 
tionally considered to be excited solutions. What can be 
seen by plotting their mass is that they play an important 
role in the presence of a gauge field and a negative cos- 
mological constant. In particular, the one-node branch in 
Fig. [3] seems to interpolate smoothly between the scalar 
dominated and gauge field dominated solutions. More 
importantly, as can be seen from Fig. ® the fact that the 
one-node and zero-node solutions come so close in terms 
of their energies seems to suggests that there is a mech- 



anism taking place at large values of the charge whereby 
the oscillatory (usually excited) solutions come very close 
to the minimal energy (ground state) solution. 
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FIG. 3: Mass as a function of scalar field at origin (left) and 
as a function of particle number (right) in the large q regime. 

A Zero Temperature Phase Transition 
The main property of a quantum phase transition [l9| is 
the existence of a value of a coupling g c at which there 
can be a level-crossing where an excited level becomes 
the ground state. Usually this creates a point of non- 
analyticity of the ground state energy as a function of 
g. Our previous graphs in Figs. [2] and [3] show this non- 
analyticity and level crossing behavior. To emphasize 
this point, we plot the critical mass of a boson star as a 
function of charge in Fig.UJ note the previously defined q\ 
are qi are the two inflexion points in this graph. Allowed 
values for boson star masses lie below and to the right of 
the curve. 

In order to interpret Fig. [4j we note that, for a given 
scalar charge q, boson stars exist for a range values of 
M. What we plot is the critical mass, namely the largest 
mass for the scalar dominated branch, and the minimum 
and the maximum masses for the gauge field dominated 
branch. To be more graphic, the right panel in Fig. [2] 
contains three distinguished points: the largest mass for 
the green branch and the smallest and largest mass for 
the red branch. These are precisely the three points in 
the phase diagram in the region q\ < q < q2- For q < qi 
there is only one type of solution and we simply plot its 
maximum mass. Similarly for q > (7 2 , we plot only the 
maximum mass. In this work we are not going to be con- 
cerned with the existence of a maximum charge, g m ax ; 
above which no regular solution exists. However, prelim- 
inary investigations suggest that such a point exists, and 
we will discuss its determination elsewhere. 

From a holographic point of view, a plot of mass versus 
the time component of the Maxwell field represents the 
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FIG. 6: The behavior of the order parameter <ti as a function 
of the charge q. 



FIG. 4: The phase diagram of charged boson stars in AdS. 



free energy as a function of chemical potential. While we 
use the gauge invariant A defined in ((6]), our asymptotic 
boundary conditions relate this to the value of u. We 
thus plot the mass versus u in Fig. [5] for solutions in the 
intermediate and large charge regimes. 
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FIG. 5: The Mass as a function of the frequency cj, which 
is a holographic proxy for chemical potential. The left panel 
corresponds to qi < q < q-2, while the right panel corresponds 
to q > q%. 



A very important role in phase transitions is played by 
the order parameter. In this case we identify it as the 
value of the scalar at asymptotic infinity, o\ . In particu- 
lar, in the case of the holographic superconductors [20j j . 
o\ corresponds to the expectation value of an operator of 
mass dimension A. In Fig. |6]wc plot o\ for critical mass 
boson stars (those following the curve of Fig. 0]), and we 
see a sharp drop around the value of q2 ■ 



Conclusions 

In the context of a charged scalar field minimally coupled 
to Einstein-Maxwell gravity with a negative cosmological 
constant, we have constructed explicit solutions and es- 
tablished a phase diagram of boson stars. The new type 
of solutions contain a gauge-field-dominated branch that 
represents a sort of "Geon" as originally envisioned by 
Wheeler, that is, a particle-like solution from mostly the 
smooth, classical fields of electromagnetism coupled to 
general relativity. 

Let us conclude by commenting on some open ques- 
tions stemming from our work. The main difference in 
gravitational collapse in asymptotically AdS spaces is the 
presence of a timelike boundary at spatial and null infini- 
ties. Under these conditions, the question of stability of 
our configurations is a particularly important one. In the 
phase space of initial conditions for gravitational collapse 
in AdS it is important to consider that a small pertur- 
bation does not escape to spatial infinity and actually 
returns to the center of AdS. We postpone such detailed 
study to a separate publication but will speculate on the 
role that boson stars might play. 

Already in the context of asymptotically Minkowski 
gravitational collapse, boson stars [21| play a particularly 
important role; they are the basis for a different type 



of critical collapse [14j. It is natural to expect that in 



asymptotically AdS spaces they will figure prominently. 
It will be interesting to pursue the construction and inves- 
tigation of properties of real boson stars, similar to those 
of asymptotically Minkowski spacetimes 21|, to asymp- 



totically AdS. One would expect that real boson stars, if 
they exist, get similarly destroyed after a certain number 
of bounces of the scalar field at spatial infinity. Recently, 
however, the conjecture that boson stars might be non- 
linearly stable has been advanced in [22| based on the 
observation that a mass scale might prevent the resonant 
turbulent mechanism of Q from being realized. 

Scalar boson stars have been previously considered in 
asymptotically AdS spacetimes by [23l - l25j . Our work has 
concentrated on charged objects but a systematic anal- 
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ysis of those configurations with the corresponding in- 
terpretation and extension to the AdS/CFT applications 
ought to be presented. We will report on such studies, 
including the neutral case, in a separate publication. 

By exploring the form of the boson star phase dia- 
gram, we have argued in favor of a gravity dual of a zero 
temperature transition for the dual field theory. It will 
be interesting, using the holographic dictionary, to com- 
pute explicitly various transport properties of such dual 
condensed matter systems. 
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